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RESUMEN 

Del 1 al 4 de febrero de 2016, se realizo en la Universidad de las Fuerzas 
Armadas ESPE, el XXII Curso Internacional de Estructuras, orientado al estudio de 
los Elementos Finitos, tema conocido a nivel mundial porque permite resolver 
problemas complejos en forma sencilla, con la ayuda de un ordenador. 

En la primera clase, Amjad Aref presento el desarrollo numerico del 
Elemento Finito Triangular, conocido como CST (Constant Strain Triangles) y 
posteriormente Enrique Morales indico el uso del programa LS-DYNA. 

En la Universidad de Fuerzas Armadas ESPE, desde el 2009, Roberto 
Aguiar viene desarrollando el sistema de computacion CEINCI-LAB para el 
analisis estatico y dinamico de estructuras, por lo que en este arttculo en primer 
lugar se presenta el marco teorico, muy bien detallado del elemento finito CST y la 
incorporacion de nuevos programas al CEINCI-LAB. Se complementa el estudio 
con la solucion de un problema estatico y otro dinamico de una viga de seccion 
constante en voladizo. 

Palabras Claves: Elemento finito CST. Tension Plana. Sistema de computacion 

CEINCI-LAB 


Recibido: Febrero de 2016 
Aprobado: Marzo de 2016 
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SUMMARY 

From February 1 st -4 th , 2016, the XXII International Course Structures was 
carried out at Universidad de las Fuerzas Armas ESPE, aimed at studying the 
Finite Element Analysis which is known worldwide because it allows to solve 
complex problems in a simple way with the aid of a computer. 

In the first class, Amjad Aref presents the numerical development of the 
Triangular Finite Element which is known as CST (Constant Strain Triangle) and 
later Enrique Morales explained the use of LS-DYNA program. 

At Universidad de las Fuerzas Armas ESPE, Roberto Aguiar has been 
developing the system of CEINCI-LAB computer for static and dynamic structural 
analysis since 2009. The first part of this paper presents a detailed theoretical 
framework of CST finite elements, and them its application through new programs 
in CEINCI-LAB. Finally, the paper presents a solution to the static and dynamic 
problem associated to a cantilever beam with a constant section. 

Keywords: finite element CST. Flat stress. System of CEINCI-LAB computer. 


1. INTRODUCCION 


El metodo del elemento finito se ha convertido en una herramienta 
poderosa en la solucion de problemas de ingenierfa. Las aplicaciones van desde 
un analisis esfuerzo y deformacion de estructuras hasta el analisis de flujo de 
calor, magnetico y otros problemas. Mediante este metodo a las regiones 
complejas se las discretiza en formas geometricas simples que son los elementos 
finitos. (Onate, 2013) 

Este artfculo es una introduccion, para quienes desean estudiar los 
elementos finitos. Por este motivo se realiza un analisis estatico y dinamico de una 
viga en voladizo, para que sea comprensible su teorfa y los programas que se van 
a desarrollar en MATLAB. 

Se destaca que para resolver una viga se puede emplear el metodo de las 
diferencias finitas (Aguiar, 2011). Es una opcion de calculo, tanto en elementos 
finitos como en diferencias la solucion de la ecuacion diferencial que gobierna el 
comportamiento de vigas a flexion se transforma en la solucion de un sistema de 
ecuaciones lineales. (Aguiar y Mroginski, 2010). 

En este artfculo se va a utilizar el elemento finito denominado: Triangulo 
de Deformacion Constante (CST), en donde la region bidimensional va a ser 
dividida en triangulos de lados rectos. Se permite que cada vertice o nodo se 
desplace en el sentido X, Y. Por lo tanto, cada nodo tiene dos grados de libertad 

(gdi). 


Para la programacion se va utilizar el software Matlab. Las coordenadas 
nodales se almacenan en un arreglo bidimensional que representa el numero total 
de nodos y sus respectivos grados de libertad, CG. Para el ensamble directo se 
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necesita un vector que contenga los grados de libertad en los nodos de cada 
triangulo, VC. Con la matriz de rigidez del elemento triangular y el VC se ensambla 
la matriz de rigidez total, Kt. A partir de la matriz de rigidez se procede a realizar 
el analisis estatico y dinamico. 


2. ELEMENTO FINITO TRIANGULAR (CST) 

La formulacion bidimensional de elementos finitos mediante el triangulo de 
deformacion constante es el mas sencillo, representa una pequena porcion de la 
estructura en 2D. Este elemento tiene los desplazamientos, las componentes de 
traccion y los valores de la fuerza de cuerpo distribuida en funcion de la posicion 
indicada por (x, y). Los nodos estan numerados en sentido anti horario, como se 
presenta en la figura 1, en que el nudo 1 es el vertice del angulo de 90 grados, al 
numerar de esta forma se garantiza que la matriz de rigidez del elementos sera 
definida positiva. 



2.1. Coordenadas reales y naturales 


Las coordenadas (x,y), del elemento finito de la figura 1, son mapeadas en 
las coordenadas uyv respectivamente, coordenadas naturales. Aref (201 6) 


u{x,y~) — a t + a 2 x + a 3 y 

(1) 

v(x,y ) = a 4 + a 5 x + a 6 y 

(2) 


Los valores de a, -a 6: de las ecuaciones 1 y 2, son parametros 
desconocidos que varfan dependiendo del elemento. Para determinar los valores 
se encuentran los desplazamientos en sentido X, en cada uno de los nudos o 
nodos de los vertices del triangulo de la figura 1 , con lo que se tiene. 
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% = u^x^y^ = a 1 + a 2 x 1 + a 3 y t 
u 3 = u(x 2 ,y 2 ) = a 1 + a 2 x 2 + a 3 y 2 
u 5 = u(x 3 ,y 3 ) = a x + a 2 x 3 + a 3 y 3 

En forma matricial, se tiene: 



1 

Xi 

yi 


W 3 = 

1 

X 2 

y 2 

r 2 

{u 5 ) 

.1 

* 3 

y 3 - 

U 3 J 


[A] 


Luego, para encontrar las constantes a t se tiene que invertir la matriz A 


{ a } 3 xl — [A]3x3{ u }3xl 


[ A ]" 1 


1 

2A e 


0 2 y 3 - * 3 y 2 ) 
(y 2 - y 3 ) 

0 3 - x 2 ) 


U 3 yi - xt y 3 ) 

(y 3 - yi) 

(x 3 - x x ) 


(xiy 2 - x 2 yj 

(yi - y 2 ) 

(X 2 - Xj) 


(3) 

(4) 


Donde, A e es el area del elemento finito triangular. 


2.2. Funciones de forma 


Las funciones de forma o de interpolacion definen la solucion de un 
problema, para un elemento de un portico piano con tres grados de libertad en 
cada uno de sus nudos, se tienen 6 funciones de forma (Aguiar, 2014). Ahora para 
el elemento triangular con los grados de libertad indicados en la figura 1 se 
tendran tres funciones de forma, las mismas que se hallan a partir de la siguiente 
ecuacion en que se supone que las deformaciones del elemento son lineales. 


u(x,y) = % + a 2 x + a 3 y 


▲ 


Sustituir con la matriz [4] 1 

u(x, y) = A^Cx, y)u x + N 2 (x,y)u 3 + N 3 (x,y)u 5 ^ 

Donde N1, N2, N3, son las funciones de forma de desplazamiento y se 
indican en la tabla 1. Se destaca que se ha remplazado los valores de a t de la 
ecuacion 3 y 4. 
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Donde A~^ son los elementos de la matriz inversa de A, que se encuentran 
en la ecuacion (4). Al reemplazar estos valores se halla en forma explicita las 
funciones de forma indicadas en la tabla 2. 


Funciones de Forma 

Ni (x,y) = (x 2 y 3 - x 3 y 2 ) + (y 2 - y 3 )x + (x 3 - x 2 )y 


N 2 (x,y) = U 3 yi - xiy 3 ) + (y 3 - yjx + (x 3 - x x )y 


N 3 (x, y) = (x x y 2 - x 2 y-d + (yi - y 2 )x + (x 2 - x ± )y 


Tabla 2 Funciones de Forma para el elemento CST 


Los valores de la funcion de forma se definen como la superficie plana, 
que se muestra en la figura 3 

En particular, Nt+IS^+Ns representa un piano con altura de 1 en los nodos 
1,2,3, es un paralelo al triangulo. 


N! + N 2 + N 3 = 1 


( 6 ) 


En este caso, se ha obtenido las funciones de forma en coordenadas 
reales. Una vez que se tienen las funciones de forma el problema practicamente 
esta resuelto. 
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N 1 ( x , y) 



N2(x,y) 3 




Figura 3 Funciones de Forma para elemento CST 


2.3. Matriz B 

Las deformaciones unitarias estan dados por: 

£= [e x ,e y ,y xy ] T 

du dN 1 dN 2 dN 3 

— — — — ~ + — — u 3 + — - — u 3 


dx dx 


dv dNi 


dx 


E y~ dy ~ ~dy~ Uz + 


dN ? 


dy 


U A + 


dx 


dNs 

dy 


du dv dN t dN 2 dN 3 dN 1 dN 2 dN 3 
Y xy , — — l~ — — — U-L + — — u 3 + — — u 5 + — — u 2 + — U4 + — — u 6 


xy dy ' dx dy 

En forma matricial, 


dy 


dy 


dx 


dx 


dx 


(£x 0,y)^ 

{6 0,y)} 3 *i = -U y (x,y)[ = 

U*y(*,y)J, 


3x1 


dN 1 

0 

dN 2 

0 

dN 3 

0 


rU 

dx 

dx 

dx 


U 2 

0 

dN 1 

0 

dN 2 

0 

dN 3 

■ 

U 3 

dy 

dy 

dy 

u 4 

dN 3 

dN 1 

dN 2 

dN 1 

dNs 

dN-L 


u 5 

- dy 

dx 

dy 

dx 

dy 

dx 


W 


{£ 0,y)} 3xl = [B(x,y)] 3 x 6 {q] 6 


( 7 ) 


( 8 ) 
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f u i\ 


(y 2 - y 3 ) 

0 

(y 3 - yi) 

0 

(yi - y 2 ) 

0 

0 

(*3 - x 2 ) 

0 

(*i - * 3 ) 

0 

(x 2 - x ± ) 

Os - * 2 ) 

(y 2 - y 3 ) 

(Xi - x 3 ) 

(y 3 - yi) 

[x 2 - xj 

(yi - y 2 ). 


u 2 

Uo 

u 4 > 
u 5 


Vu 6 J 


( 9 ) 

De donde la matriz de compatibilidad de deformaciones que sirve para 
calcular la matriz de rigidez es la siguiente. 


y 23 

0 

y 3i 

0 

Vl2 

0 ■ 

0 

X 32 

0 

X 13 

0 

X 21 

- x 22 

y 23 

x 13 

y 3 1 

X 21 

yu 


( 10 ) 


Notese que la matriz B esta completamente definida, se puede obtener 
directamente para un elemento finito cualquiera, esto facilita notablemente el 
calculo de la matriz de rigidez del elemento. 

2.4. Matriz de rigidez del elemento CST 

La matriz de rigidez del elemento en forma general viene definida de la 
siguiente manera: 


K = 


" 


B T D Bdv e — t 


" 


B T D BdA 


( 11 ) 


Normalmente la integracion se realiza utilizando la cuadratura de Gauss 
(Aguiar 2016, Aref 2016) pero para el caso del elemento finito CST, la matriz de 
rigidez del elemento se halla con la siguiente expresion. 

K= A e t[B] T [D][B] 


Donde: K es la matriz de rigidez del elemento finito CST; A e es el area del 
elemento finito; t es el espesor del elemento; B es la matriz de compatibilidad 
presenta en la ecuacion (10); D es la matriz que relaciona el vector de esfuerzos 
con el vector de deformaciones. 
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Para el caso piano esta matriz es de tres por tres y para el caso de 

tension plana, se tiene: 

ri v o 


D = 


1 — V 2 


v 1 0 

1—17 

0 0 


2 


(13) 


Donde E es el modulo de elasticidad del material; v es el modulo de 


Poisson. 


3. PROGRAMAS PARA EL ANALISIS ESTATICO 

Para la utilizacion del programa, la numeracion de los nudos, elementos y 
grados de libertad se indican en la figura 3. 



Figura 3 Grados de libertad, Numeracion de nudos y triangulos. 


Notese que primero se han numerado los nudos restringidos, luego los 
restantes nudos de la malla, de abajo hacia arriba. Los elementos se han 
numerado de izquierda a derecha y de abajo hacia arriba. Finalmente se indica 
que cada nudo tiene dos grados de libertad un desplazamiento horizontal y un 
desplazamiento vertical; como se va a resolver un problema estatico los grados de 
libertad se numeran de corrido. 

La matriz de rigidez de la estructura, se halla por ensamblaje directo y 
para ello se necesita determinar la matriz de rigidez de cada elemento finito y los 
vectores de colocacion, que no son nada mas que los grados de libertad de los 
nudos, escritos en el orden en que fue deducida la matriz de rigidez del elemento. 

En la figura 4 se aprecia con dfgitos mas grandes la numeracion de los 
nudos y con dfgitos mas pequenos la numeracion de los nudos del elemento finito, 
teniendo en cuenta que se numeran en forma anti horaria, empezando en el 
vertice de angulo recto. 
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Figura 4 Numeration de nudos en sentido anti horario para cada elemento. 

Con estas acotaciones a continuacion se presentan los nuevos programas 
que han sido incorporados al sistema de computacion CEINCI-LAB para el 
analisis con elementos finitos CST. 

■ dibujo_elemFinito_triangulo( x ,y , L , h ) 

Esta funcion sirve para graficar todas las divisiones triangulares. 

x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

L Longitud total, en metros, 

h Altura del elemento. 


■ [K]=Kmiembro_ElemFinito_triangulo (L,e,h,x,y,u,E) 

K Matriz de rigidez del elemento triangular. 

E Modulo de elasticidad, en T/m2. 

u Poisson, 

e Espesor. 

h Altura del elemento. 

x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

■ [C G] = cg_elemFinito_ triangulo(x,y) 

CG Matriz de coordenadas generalizas. 

x Division en el eje x al elemento. 

y Division en el eje y al elemento. 


■ [VC]= vc_elemFinito_triangulo(x, y,CG) 

VC Matriz que contiene vector de colocacion. 
x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

■ [K[]=krigidez_elemFinito_triangulo{x,y,\/C,K) 

Kt Matriz de rigidez de la estructura. 
x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

VC Matriz que contiene vector de colocacion. 

K Matriz de rigidez del elemento triangular. 

■ desplazamien fo(x,y,L,h,q) 

Esta funcion sirve para graficar los desplazamientos de los grados de 
libertad verticales. 
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x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

L Longitud total, en metros, 

h Altura del elemento. 

q Matriz que contiene los desplazamientos. 


4. PROGRAMA PARA EL ANALISIS DINAMICO 


Para el analisis dinamico es necesario numerar al primero o al ultimo las 
coordenadas principales. En este caso se numeran al comienzo, como se indica 
en la figura 5. Se esta considerando que en sentido vertical se tiene un grado de 
libertad. Con esta numeracion no se puede aplicar la triangularizacion de Gauss 
para hallar la matriz de rigidez condensada a las coordenadas principales. 



Figura 5 Grados de libertad, Numeracion de nudos y triangulos. 


Ahora se debe encontrar la matriz que contiene a los vectores de 
colocacion de cada uno de los elementos, para ello se mantiene la misma notacion 
del problema estatico pero para que sea didactico se repite la figura 4 en la figura 
6 . 



Figura 6 Numeracion de nudos en sentido anti horario para cada elemento. 

Los programas, en su mayor parte son los mismos que para el analisis 
estatico sin embargo se los presenta nuevamente, para que sean agrupados en 
dos carpetas, la una para analisis estatico y la otra para analisis dinamico. 

■ dibujo_elemFinito_ triangulo{ x ,y , L , h ) 

Esta funcion sirve para graficar todas las divisiones triangulares. 

x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

L Longitud total, en metros, 

h Altura del elemento. 
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■ [K]=Kmiembro_ElemFinito_triangulo (L,e,h,x,y,u,E) 

K Matriz de rigidez del elemento triangular. 

E Modulo de elasticidad, en T/m2. 

u Poisson, 

e Espesor. 

h Altura del elemento. 

x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

■ [C G] = cg_elemFinito_ triangulo_sismo(x , y ) 

CG Matriz de coordenadas generalizas. 

x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

■ [VC]= vc_elemFinito_triangulo{x, y,CG) 

VC Matriz que contiene vector de colocacion. 
x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

■ [Kt gti\]=krigidez_elemFinito_triangulo_sismo(x,yyC,K) 

Kt Matriz de rigidez 

x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

VC Matriz que contiene vector de colocacion. 

K Matriz de rigidez del elemento triangular, 
gdl Numero de grados de libertad. 

■ [T,fi,OM]=orcte/7_e/g(KL,M) 

KL Matriz de rigidez lateral 

M Matriz de masas. 

T Periodo de vibracion. 

Fi Modos de vibracion. 

OM Frecuencia 


5. APLICACIONES 

5.1. Viga en voladizo 

Se propone calcular el desplazamiento que genera una fuerza de 5T, 
aplicada en el extremo como se indica en la figura 7, a una viga en voladizo de 
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longitud de 3 metros. Mediante el uso del programa de elemento finito triangular y 
comparar el resultado mediante la teorfa de estructuras. 

Seccion: Base=40cm Material: E=240000kg/cm 2 

Altura=30cm 




h=0.40m 



A continuacion se indica el programa que se debe elaborar para resolver el 
problema indicado en la figura 7. Utilizando el sistema de computacion CEINCI- 
LAB. 


IDatos 

L=3; llongitud en metro 

h=0.4; %Altura en metro 

e=0.3; %espe3or en metro 

x =40; lx division de los nudos verticales 

y =10; ly division de los nudos horizontales 

E=240Q000; It/m2 

u=0.2; Ipoisson 

I 

dibujo_elemFinito_triangulo (x,y,I,h) ; 
[K]=Kmiembro_ElemFinito_triangulo (L,e,h,x,y,u,E) ; 
[CG]=cg_elemFinito_triangulo (x,y) ; 

[VC] =vc_elemFinito_triangulo (x, y, CG) ; 
[Kt]=lcrigidez_elemFinito_triangulo (x,y,VC,K) ; 
gdl=2* (x+1) * (y+2) ; 

P=zeros ( gdl , 1 ) ; 

P(gdl)=-5; 

q=Kt\P; 

desplazamiento(x,y,L,h,q) ; 

I Calculo por teoria de estructuras 

I=0.3*0.4 A 3/12;EI=E*I;L=3.0;P=5; 

des=P*L A 3/(3*EI); 

Figura 8 Lfneas de Programacion. 

En la figura 9 se indica el mallado con el cual se ha resuelto el problema y 
en la figura 10 se tiene los desplazamientos verticales. Finalmente, se compara el 
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desplazamiento vertical en el voladizo con el que se obtiene con la teoria de 
estructuras, se observa que existe una muy buena aproximacion. 



Figura 9 Mallado de triangulos. 



Figura 10 Desplazamiento en la ordenada 

Resultados 

Elementos finitos: q = 0.01 1 2 m 
Teona de estructuras: q = 0.01 1 7 m 


5.2. Propiedades dinamicas de viga en voladizo 

Se propone calcular los periodos de vibracion de la viga en voladizo, 
presentada en la figura 11 que tiene una longitud de 3 metros. Mediante el uso del 
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programa de elemento finito. Comparar el resultado utilizando el elemento finito Q4 
y mediante la solucion de una viga modelada como un sistema continuo. La carga 
que gravita de 2 T/m sirve para determinar la matriz de masas 


Seccion: 


Base=40cm Material: E=240000kg/cm 2 

Altura=30cm 



Figura 11 Viga en voladizo. 

Las coordenadas principales son los desplazamientos verticales y se halla 
la matriz de Masas, asociada a esos grados de libertad obteniendo la masa total y 
dividiendo para el numero de coordenadas principales. A continuacion se presenta 
el programa desarrollado utilizando los programas de CEINCI-LAB. 
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%Datos 

L=3; %longitud en metro 

h=0 . 4 ; %Altura en metro 
e=0.3; %espesor en metro 

x=200; %x division de los nudos verticales 

y=5; %y division de I 03 nudos horizontales 

£=2400000; %t/m2 
u=0.2; %poisson 

% 

dibujo_elemFinito_triangulo (x, y,L,h) ; 

[K] =Kmiembro_ElemFinito_triangulo (L,e,h,x,y,u,E) ; 

[CG] =cg_elemFinito_triangulo_sismo (x, y) 
[VC)=vc_elemFinito_triangulo (x,y,CG) ; 

[Kt gdl] =krigidez_elemFinito_triangulo_sismo (x, y, VC, K) ; 

% K lateral 

na=x+l ; ngl=gdl ; 
kaa=Kt (1 :na, 1 :na) ; 
kab=Kt ( 1 : na , na+1 : ngl ) ; 
kba=kab ' ; 

kbb=Kt (na+l:ngl,na+l:ngl) ; 

KL=kaa-lcab* inv ( kbb) *kba; 

% Anali3is dinamico 

MT=2*L; m=WT/( (x+1) *9.8) ; 

M=zeros (l,x+l) ; 

□ |for i=l : x+1 

M(l,i)=m; 

L end 

M^diag (M) 

[T, fi,OM]=orden_eig (KL,M) 

Figura 12 Linea de Programacion. 

Para tener una buena aproximacion en el calculo de los perfodos de 
vibracion es necesario trabajar con muchos elementos finitos, lo que no sucede 
cuando se trabaja con el elemento finito Q4. En la figura 14 se presenta la malla 
en que se ha considerado 100 divisiones en sentido X, y 5 divisiones en sentido Y. 

Se obtuvo la matriz de rigidez de rigidez asociada a todos los grados de 
libertad y se condenso a las coordenadas verticales, como se observa en la figura 
12. Los primero cinco perfodos de vibracion se presentan en la tabla 3. 
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Figura 13 Mallado de triangulos. 


Modos 

1 

2 

3 

4 

5 

Periodos de 
Vibracion 
(segundos) 

0.1134 

0.0192 

0.0074 

0.0042 

0.0028 


Tabla 3 Periodos de Vibracion de los cinco primeros valores. 


• Elemento Finito Q4 

En la figura 14 se presenta la discretizacion del dominio al trabajar con el 
elemento finito Q4 suavizado por flexion. Notese que se tiene 5 divisiones en X, y 
3 divisiones en Y. Se ha numerado solo las coordenadas verticales, que se 
considera igual a todo lo largo de la Ifnea vertical. (Aguiar, 2016) 



Figura 14 Coordenadas Principales 

Las coordenadas secundarias son los desplazamientos en cada uno de los 
puntos de la malla, de tal manera que se tienen 25 grados de libertad, se halla la 
matriz de rigidez y se condensa a 5 por 5. La matriz de masas es de5 por 5. 

La numeracion de los elementos se muestra en la figura 15 y en la tabla 4 
se indican los periodos de vibracion. 
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Modos 

1 

2 

3 

4 

5 

Periodo de 
Vibracion 
(segundos) 

0.1385 

0.0219 

0.0078 

0.0042 

0.0028 


Tabla 4 Periodo de Vibracion mediante Q4 


El periodo de vibracion encontrado con el elemento finito Q4 es 0.1385 
segundos y el que se encuentra con el elemento triangular CST es 0.1 134 seg. Se 
trabajo con pocas divisiones con el elemento finito Q4 para entender mejor la 
solucion del problema. 

• Sistema Continuo 

La viga en voladizo de flexion se puede modelar mediante un sistema 
continuo descrito con detenimiento en Aguiar (2012). Ahf se aprecia que la 
solucion de la ecuacion diferencial en derivadas parciales conduce a la siguiente 
ecuacion. 


1 + cos p Cosh p — 0 


(14) 


p = a L 

, EI 
Wn 2 — — a' 
m 


T = 


2 n 

w n 


(15) 

(16) 


(17) 
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Donde L es la longitud de la viga en voladizo; E es el modulo de 
elasticidad del material; / es el momento de inercia de la viga; in es la masa por 
unidad de longitud. Las cinco primeras raices de la ecuacion (14) son: 


Pi=1 .875; p 2 =4.694; p 3 =7.854; p 4 =1 0.996; p 5 =14.137 (18) 


Para calcular el periodo fundamental asociado al primer modo se tiene: 

Pi = 1.875; a = 0.625; W n = 53.5826; T = 0.1173 seg 

En la tabla 5 se indican los periodos para los cinco primeros modos de 
vibracion. 


Tabla 5 Periodo de Vibracion mediante un sistema continuo. 


Modo 

1 

2 

3 

4 

5 

Periodo de 
Vibracion 
(segundos) 

0.1173 

0.0187 

0.0067 

0.0034 

0.0021 


El periodo fundamental de vibracion hallado con el elemento triangular 
CST es 0.1134 seg., que es muy parecido a la solucion exacta hallada con el 
sistema continuo de 0.1 1 73 seg. 


6. CONCLUSIONES 

Se ha presentado el marco teorico que conduce al calculo de la matriz de 
rigidez de un elemento finito triangular CST, siguiendo el desarrollo numerico 
presentado en Aref (2016). Por otra parte se han desarrollado varios programas de 
computacion para el analisis estatico y dinamico de vigas en voladizo con 
elementos triangulares y se resolvieron dos ejemplos. Del estudio realizado, se 
desprende: 

■ La contribucion al sistema de computacion CEINCI-LAB para la modelacion 
de una viga en voladizo mediante el triangulo de deformacion unitaria 
constante (CST) dan resultados confiables para manipular pero se tiene que 
tener mucho cuidado debido que la malla que se va a generar debe ser muy 
densa para que estos resultados sean buenos, este principio se aplica al 
analisis estatico como al dinamico. Al realizar la modelacion con una malla 
poco densa los resultados son distorsionados y por tanto no se debe utilizar 
para un posterior calculo. 
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